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Abstract 

Using first-principles calculations of graphene having high-symmetry distortion or de- 
fects, we investigate band gap opening by chiral symmetry breaking, or intervalley mixing, 
in graphene and show an intuitive picture of understanding the gap opening in terms of local 
bonding and antibonding hybridizations. We identify that the gap opening by chiral symme- 
try breaking in honeycomb lattices is an ideal two-dimensional (2D) extension of the Peierls 
metal-insulator transition in ID linear lattices. We show that the spontaneous Kekule distor- 
tion, a 2D version of the Peierls distortion, takes place in biaxially strained graphene, leading 
to structural failure. We also show that the gap opening in graphene antidots and armchair 
nanoribbons, which has been attributed usually to quantum confinement effects, can be under- 
stood with the chiral symmetry breaking. 

*To whom correspondence should be addressed 
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High carrier mobility makes graphene promising as a channel material for field effect tran- 
sistors. 1 ' 2 The absence of a band gap in graphene, however, gives rise to poor on-off ratios in 
the transistor performance. 3 Several approaches to open a band gap have been suggested, such 
as interactions with substrates, 4 patterning into nanoribbons, 5 " 8 quantum dots 9 or into periodic 
structures containing carbon vacancies, called antidots,— ~— periodic modulation with hydrogen 
adsorption,— or imposing strain.-^— The origin of gap opening in those structures has been usu- 
ally attributed to sublattice symmetry breaking, 18 magnetic effects, 7 or quantum confinement ef- 
fects.- 59 -""— 1 ^ 9 However, it is unclear how quantum confinement effects open a band gap in graphene. 
In fact, there are many counterexamples. Under the time reversal symmetry, i.e., without including 
magnetic effects, graphene nanoribbons with zigzag edges do not have a band gap despite ID con- 
finement. 5 " 8 Even with armchair edges graphene nanoribbons of certain widths have a zero band 
gap when lattice distortions at edges are disregarded. 5-7 Further, studies on graphene under peri- 
odic potentials found that confinement itself cannot open a band gap.— ^ A clear understanding 
of the gap opening is necessary to accelerate achieving practically viable gapped graphene. 

Graphene has a zero band gap with four half-filled degenerate states at the intrinsic Fermi 
level. 22 The fourfold degeneracy, consisting of two degenerate states at two nonequivalent Dirac 
points (K and K'), comes from the crystal symmetry of graphene's honeycomb lattice. Band gap 
opening in graphene thus implies breaking of the symmetry. Analytical studies based on effec- 
tive Hamiltonians-2 3 - - — have proposed various symmetry-breaking mechanisms of gap opening, in- 
cluding sublattice symmetry breaking, chiral symmetry breaking, spin-orbit coupling, etc. Among 
those, chiral symmetry breaking, or intervalley mixing, that couples Bloch states at the two Dirac 
valleys with each other has fundamental importance. Since the electronic states at the K and K' 
valleys in graphene represent the massless Dirac fermion spectrum of different chiralities 22 (with 
the spin being the pseudospin defined in the sublattice space), the gap opening by chiral symme- 
try breaking corresponds to the mass gap generation of the massless Dirac fermions in quantum 
electrodynamics (QED). 28 

In this paper, we study the gap opening by chiral symmetry breaking in graphene using first- 
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principles calculations. We examine in detail the electronic wavefunctions of gapped graphene 
having high-symmetry distortion or defects, and show that gap opening by chiral symmetry break- 
ing can be understood easily in terms of local bonding and antibonding hybridizations. Especially 
we identify that the chiral symmetry breaking in honeycomb lattices via electron-lattice coupling 
is an ideal 2D manifestation of the ID Peierls metal-insulator transition and show that spontaneous 
2D lattice distortion occurs in graphene when biaxial strain is applied, initiating structural fail- 
ure. Finally we show that the gap opening in graphene antidots and armchair nanoribbons can be 
understood due to chiral symmetry breaking. 



Results and Discussion 

We begin by discussing how gap opening occurs in graphene within time reversal symmetry. In a 
nearest-neighbor tight-binding model, 22 the electronic structure of pristine graphene is described 
by two bands, each from the bonding (jz) and antibonding {%*) hybridization of two sublattice 
atoms' p z orbitals ( [figure] [l][]lfr ). The two bands become degenerate at K and K'. The four de- 



generate states, 7z(K), 7Z*(K), 7z(K'), and 7Z*(K'), can undergo energy splitting in two ways, 26 the 
7Z-7Z* and K-K' mixing. The iz-iz* mixing or sublattice symmetry breaking occurs by asymmet- 
ric on-site energies for the two sublattices, producing dehybridization into sublattice Bloch states, 
e. g. , at the K point, y/^ {K + q) = ~y n {K + q) ± yA^* (K + q). The K-K' mixing or chiral symmetry 
breaking occurs by periodic potentials that provide a momentum transfer of Q = ±2K (equivalently 
Q = ^K). Among all four K-K' mixing possibilities of 7Z(K)-7z(K'), 7Z*(K)~7Z*(K'), 7Z(K)-7Z*(K') 
and 7Z*(K)-7z(K'), the first two play a role in gap opening, because only those mixings introduce 
level crossing at the intrinsic Fermi level, as depicted in |[figure][l][]lb for the tz(K)-tz(K') mix- 



ing. In an effective 4x4 Hamiltonian acting on a four-component spinor consisting of the four 
degenerate states, those gap-opening mixing terms correspond to the mass terms in the Dirac equa- 
tion. They become finite in graphene's honeycomb lattice upon the Kekule distortion. 23 > 24 
The Kekule-distorted structure ( [figure][l][]lj :) has a Wigner-Seitz cell containing a single ben- 



zene unit (the gray hexagon in [figure] [1 ][]!{ :), where the C-C bonds consist of alternating short 
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Figure 1: Band gap opening by chiral symmetry breaking in graphene's honeycomb lattice, (a) 
Band structure of pristine graphene along the line connecting the Brillouin zone center (T) and 
the K point. The n{K)-n*{K) mixing can occur by sublattice- symmetry breaking potentials. The 
upper and lower panels show wavefunctions at T, K and K'. The wavefunction coefficients of C 2p z 
orbitals are reflected in the radius and color of a circle at each atomic site, with the phase angle of 
0, 2n/3, and 4n/3 being represented by red, green, and blue, and the intermediate angles by their 
interpolation, (b) n(K)-n(K') mixing or chiral symmetry breaking. A level crossing occurs by the 
momentum transfer of Q = ±2Kx by perturbing potentials, (c) Atomic structure of graphene with 
the Kekule distortion. The thick bonds are shorter than the thin bonds. Two equivalent Wigner- 
Seitz cells are shown, (d) Calculated band structure of a Kekule-distorted structure with ~ 7 % 
bond length asymmetry. (e,f) Wavefunctions of the lowest empty states. (g,h) Wavefunctions of 
the highest occupied states. 
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and long bonds. It has been studied earlier as a possible low-energy structure for aromatic hy- 
drocarbons 29 and carbon nanotubes,— i2ii22a2J. involving the Peierls instability— as in ID atomic 
chain, but such an instability was found not to occur in 2D graphene at and far below room temper- 
ature.—^ (We note that the geometry depicted in [figure][l][]l{ ; is actually the inverse distortion 
of the typical Kekule distortion in earlier studies,— for which we will discuss shortly.) 
We performed first-principles calculations for a forced Kekule-distorted structure to scruti- 



nize chiral symmetry breaking in graphene. The calculated band structure ([figure][l][]l i) shows 
that the distortion breaks the fourfold degeneracy of pristine graphene into two twofold degenera- 
cies, opening a band gap. The wavef unctions of the highest occupied and lowest empty states 
( [figure][l][]lfc -h) reflect the K-K' mixing: They are the mixed states of the four degenerate 
states of pristine graphene, \]/7i gu (q) = Y7t(K + q) ± \j/n{K' + q) by the n(K)-n(K') mixing and 
\\r % * g {q) = \ir n *(K + q) ± yf K *(K' + q) by the n* {K)-%* (K') mixing. The wavefunctions of the 
chirality-mixed states are the standing waves made of plane waves moving along K and K' (= — K) , 
and show a characteristic feature that every C-C pair has a phase difference of and n, exhibiting 
apparent local bonding or antibonding hybridization. 

The calculated electronic wavefunctions show two interesting features of chiral symmetry 
breaking in honeycomb lattices. One is that the Kekule distortion in honeycomb lattices is an 
ideal 2D extension of the ID Peierls distortion. When we take an alternative Wigner-Seitz cell (the 



red hexagon in [figure] [1][]1 the distortion shortens all the intracell bonds, while it lengthens 
all the intercell bonds, as the Peierls distortion does in a two-atom unit cell atomic chain. 33 Fur- 
ther, the gap opening occurs in the same manner in terms of the modification of local bonding and 
antibonding hybridizations. Among the four chirality-mixed states {7t g>u , Ki,g} that are degenerate 
in pristine honeycomb lattices, the n u and Tt* states get a lower energy by the distortion because 
it enhances bonding hybridizations in the intracell bonds while it lessens antibonding ones in the 
intercell bonds ( |[figure][l][]lg ,h). For the 7t g and n* states, on the other hand, the opposite occurs 
( [figure][l][]lfc ,f), raising their energy. We found that there is more correspondence between the 
2D honeycomb and ID linear lattices: Both lattices in nearest-neighbor tight-binding models show 
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Figure 2: Spontaneous chiral symmetry breaking in graphene under biaxial strain, (a) The energy 
profile calculated as a function of the relative intracell atomic distances with the equibiaxial strain 
e = £ x = £ y of 14, 15, and 16%. The upper pannels show the schematic atomic geometries of 
two local minima, A and B, corresponding to the inward and outward Kelule distortions, (b) The 
intracell and intercell bond lengths for the inward distortion (A), as a function of £. They bifurcate 
at e = 14.2%. ao is the equilibrium bond length at e = 0. (c) The resulting band gap for the inward 
distortion (A). 

massless Dirac fermion spectra with fourfold and twofold degeneracy, respectively, at the intrinsic 
Fermi level and both undergo two common gap opening mechanisms originating from asymmetry 
in on-site energies and that in hopping integrals, respectively, each corresponding to sublattice and 
chiral symmetry breaking (Supporting Information). This indicates that, regarding the electronic 
structures within tight-binding descriptions, the honeycomb lattice, not a square or rectangular lat- 
tice, is the 2D extension of a linear lattice, and to the same extent the Kekule distortion is the 2D 
extension of the Peierls distortion. 

Another interesting feature is the parity symmetry revealed in the chirality-mixed states. Ac- 
cording to QED, massive Dirac fermions have an intrinsic parity, pertaining to internal structures 
of point particles, similar to the intrinsic spin, and their antifermions have an opposite parity. 34 The 
wavefunctions in [figure][l][]lfc -h have parity symmetry and show opposite parities for the occu- 



pied (n u , 7t*) and empty (7Z g ,7C*) states. The parity symmetry is defined within the Wigner-Seitz 
cell having six atoms inside, and each of the Wigner-Seitz cells that constitute their own hexagonal 
lattice has a constant wavefunction amplitude throughout the 2D space. The wavefunctions there- 
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Figure 3: Tearing of a graphene sheet via spontaneous Kekule distortion under biaxial strain. 
First-principles molecular-dynamics calculations are performed for graphene under an equibiaxial 
strain e = 15.5% at T = 300 K. A (6 x 6) supercell is used (red rhombus in the figures), (a) t = 20 
fs. With the applied strain, graphene undergoes a prompt Kekule distortion, (b) t = 50 fs. By 
thermal fluctuation, some longer bonds appear and trigger weakening of neighboring long bonds, 
(c) t = 100 fs. Boundaries made of longer bonds become prominent, (d) t = 300 fs. Each domain 
of graphene contracts to relieve the tensile stress and recover equilibrium atomic distances. 



fore visualize the product of the parity wavefunction and the spatial wavefunction of massive Dirac 
fermions at rest in a hexagonal lattice. The two degenerate states for the empty (occupied) states 
correspond to different intrinsic spin states of the massive Dirac fermions (antifermions); The spin 
states are seen as the relative phases between sublattice atoms in the wavefunctions (i.e., the pseu- 
dospin for graphene). Thus the present results show that the internal structure of 2D point particles 
in a continuum description is visualized in a lattice description as a six-site internal structure of a 
hexagonal lattice. 

The band gap opening by the symmetry-breaking distortion that normally does not occur in 2D 
graphene takes place spontaneously in biaxially strained graphene. Our first-principles calcula- 
tions show that a spontaneous Kekule distortion occurs at equibiaxial strain of £ = e x = £ y > 14% 



([figure] [2] []2i,b). It shows asymmetric two local minima, the inward (A) and outward (B) dis- 
tortions. The outward Kekule distortion leading to dimerization of carbon atoms, which has been 
discussed previously 23 ' 24 ' 29 ~ — as a 2D extension of the Peierls distortion, is less stable than the 
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Figure 4: Band gap opening in graphene by hexagonal antidot superlattices. (a) Atomic structure 
for the (4 x 4) supercell (N = 4). The antidot consists of six carbon vacancies in the shape of 
benzene ring. Dangling bonds are passivated by hydrogen atoms so that the sp 2 bond network 
is intact, (b) Calculated band gaps as a function of the antidot supercell size (AO- Red and blue 
dots represent different characteristic of the data. Solid lines are fitted to N~ 2 , proportional to the 
density of defects. 

inward distortion. The preference to the inward distortion of graphene comes from the o bonds in 
graphene: We performed calculations of a comparative honeycomb lattice made of hydrogen atoms 
and found that it favors the outward distortion (Supporting Information). The calculations also in- 
dicate that the strong a bonds are responsible for the lack of the Kekule distortion in graphene 
at zero strain. The distortion under biaxial strain produces a finite band gap ( [figure] [2] []2{ :) that 



grows rapidly with increasing strain after the onset at e = 14.2% up to £ = 17.8% where the a* 
band starts to descend below the empty % band. This chiral-symmetry-broken state of graphene, 
however, is subject to a structural failure triggered off at fluctuating long bonds: Our molecular 
dynamics simulations of graphene under a biaxial strain of 15.5% result in tearing of the graphene 



sheet via the Kekule distortion ( [figure] [3] []3 ). We note that the symmetry breaking and gap open- 
ing under biaxial strain is distinct from that by uniaxial strain,-^^ where the gap opening occurs 
by the gradual merging of the two Dirac points caused by the uniaxial distortion of the honeycomb 
lattice.— 

The chiral symmetry breaking in graphene can occur not only by spontaneous lattice distortions 
but also by explicit lattice defects. We now show the manifestation of chiral symmetry breaking in 
previously reported structures of gapped graphene, for which gap opening was attributed to quan- 
tum confinement effects or others. We consider the periodic antidot defects first.— - — Hexagonal 
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Figure 5: Explicit chiral symmetry breaking in graphene by hexagonal antidot superlattices. (a) 
Atomic structure for the (6 x 6) supercell. (b) Calculated band structure. For the (6 x 6) supercell, 
the two Dirac points fold into the T point in the reduced Brillouin zone. The fourfold degeneracy 
is broken into two twofold degeneracies. (c,d) Wavefunctions of the lowest empty states, each 
corresponding to n g and n*. (e,f) Wavefunctions of the highest occupied states, each corresponding 
to K u and 71*. 



antidot superlattices with antidots composed of six carbon vacancies ( [figure] [4] []4fr ) reduce the 
translational symmetry, yet preserving graphene's C§ v point group symmetry. The calculated band 
gaps for the (N x N) supercell ( [figure] [4] []4b ) show that they are nonzero and proportional to 
the density of defects only when the size of the antidot superlattice, N, is a multiple of 3. This 
supercell-size dependence invalidates quantum confinement effects as the origin of gap opening. 



The calculated wavefunctions ( [figure] [5] [] 5 1 show that chiral symmetry breaking instead is the 
origin: The lowest empty and highest occupied states represent nothing but the chirality-mixed 
states, {iCg, K*} and {tc u , tc*}, respectively. The gap opening, or the degeneracy lifting of the four 
chirality-mixed states, can be understood as the result of a systematic removal of bonding neighbors 
for {n g , it*} and of antibonding neighbors for {n u , n*} by the antidot formation. This gap opening 
depends on the commensurability of the chirality-mixed states with the defect lattice: When they 
are incommensurate, i.e., when N is not a multiple of 3, the antidot formation removes both bond- 
ing and antibonding neighbors, resulting in a net energy change of zero, for all the chirality-mixed 
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Figure 6: Band gap opening by chiral symmetry breaking in armchair graphene nanoribbons. 
(a) Atomic structure of an armchair nanoribbon with nine carbon dimers along the width, (b) 
Calculated band structure. The fourfold degeneracy is fully lifted. (c,d) Wavefunctions of the 
lowest and second lowest empty states at T, corresponding to Tt g and it*. (e,f) Wavefunctions of the 
highest and second highest occupied states at T, corresponding to n* and % u . The edge truncation, 
depicted in the upper panel, makes the n* and % u states deviate significantly from zero energy 
and thus mix strongly with high-energy states. The edge dimer length is shortened spontaneously 
by 3.5 % compared to that in bulk. 7 It is the consequence of the stronger bonding hybridization 
between the edge dimer atoms of the second highest occupied state, compared to the antibonding 
hybridization of the highest occupied state. The edge dimer distortion thus reduces the band gap 
slightly. 



states. 

Similarly one can understand the gap opening in quasi- ID graphene nanoribbons with arm- 
chair edges,- - - where the armchair edges impose the chiral symmetry breaking. As shown in 
[figure] [6][]6| the edge truncation removes bonding neighbors for K g and n*, while it removes 



antibonding ones for n* and n u , and thereby lifts the degeneracy among them. For the K g and n* 
states ( [figure] [6][]6{ :,e) that have weak hybridizations at the edges, the energy change is small, but 



for the n* and % u states ( [figure] [6] []6{ i,f) that have strong hybridizations at the edges, the energy 



change is large, leading to strong mixing with high-energy states. Similar electronic structures can 
be produced in 2D graphene with periodic line patterns of hydrogen adsorption along the armchair 
direction (Supporting Information), where the passivation of p z orbitals by hydrogen adsorption 
plays the role of edges. The inverse proportionality of the band gap to the ribbon width, which 
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Figure 7: Band gap opening by chiral symmetry breaking in armchair graphene nanoribbons. 
(a) Atomic structure of an armchair nanoribbon with eight carbon dimers along the width, (b) 
Calculated band structure. The fourfold degeneracy is fully lifted. (c,d) Wavefunctions of the 
lowest and second lowest empty states at T, corresponding to n* and lt g . (e,f) Wavefunctions of the 
highest and second highest occupied states at T, corresponding to % u and it*. The edge truncation, 
depicted in the upper panel, makes the 7t g and it* states deviate significantly from zero energy 
and thus mix strongly with high-energy states. Here the it g state gets a higher energy despite the 
removal of antibonding neighbors for the edge atoms. It is because that the edge truncation also 
removes strong bonding pairs that make a higher contribution to the energy. Similarly the it* 
state gets a lower energy. For the K* and it u states the edge truncation removes just nonbonding 
neighbors and their degeneracy is lifted by spontaneous distortion of edge dimers. 

has been regarded as an evidence for the manifestation of quantum confinement effects for arm- 
chair nanoribbons,- 1 ^ can thus be understood as just the result of the decreasing density of chiral 
symmetry breaking defects as the width increases. 

An interesting case arises when the armchair nanoribbons have a certain width of (3n — 1) 



carbon dimers; the n = 3 case is shown in [figure] [7] []7 With these widths, the Tt g and it* states 



( |[figure][7][]7| i,f) undergo large energy splitting by the edge truncation, but the it* and lt u states 
( [[figure] [7] []7b ,e) are yet degenerate at zero energy, 5 ' 6 because they have nodal lines along the 
armchair edges and the edge truncation removes just nonbonding neighbors at both edges. The 
degeneracy is lifted actually by spontaneous dimer distortion at edges leading to 3.5 % reduction 
in the bond length. 7 This so-called edge effect is again the manifestation of chiral symmetry break- 
ing: The edge distortion enhances bonding hybridizations for it u , whereas it enhances antibonding 
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hybridizations for it* leading to the energy splitting between them. This spontaneous gap open- 
ing, allowed in graphene for the reduced elastic energy cost from the low symmetry at edges, is 
identical to the ID Peierls transition. It is especially true for the n = 1 case which is just the ID 
atomic chain system of czs-polyacetylene, 35 with only the n* and % u states being involved for the 
gap opening. This demonstrates that the gap opening by chiral symmetry breaking in graphene is 
indeed a 2D superset of the ID Peierls transition. 

Conclusion 

We have shown using detailed analysis of the electronic wavefunctions of gapped graphene that the 
gap opening by chiral symmetry breaking in graphene can be understood as an ideal 2D superset 
of the ID Peierls transition and also as the degeneracy lifting of four-chirality-mixed states, which 
are represented as the network of local bonding and antibonding hybridizations. Our study have 
shown that this understanding is useful, providing a simple, unified description of the gap opening 
in 2D graphene antidots and quasi- ID armchair nanoribbons and predicting structural failure of 
biaxially strained graphene via spontaneous 2D lattice distortion. 

Methods 

Our first-principles calculations are based on the density functional theory (DFT) employing the 
generalized gradient approximation and the projector-augmented- wave method as implemented in 
VASP.— «2Z Valence electronic wavefunctions are expanded in a planewave basis set with a cutoff 
energy of 280 eV. In our supercell calculations, graphene layers are separated from each other by 
8 A. The fc-point integration was made at a uniform fc-point mesh of (30 x 30) per unit cell. The 
atomic positions are relaxed until residual forces are less than 0.02 eV/A. 
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